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Using mainly analytical arguments, we derive the exact relation r^max = y ^ for the max- 

imal vacuum value of the Higgs field for static gravitational global monopoles. For this 
value, the global monopole bifurcates with the de Sitter solution obtained for vanishing 
Higgs field. In addition, we analyze the stability properties of the solutions. 
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Figure 1: Schematic of the space of static solutions. For large ff, only de Sitter exists 
and is stable. As ff decreases, a bifurcation occurs and de Sitter exchanges stability 
with the static monopoles. 



In a recent publication one of us (SL) presented numerical results on static, 
self-gravitating global monopoles. Considered as topological defects, global monopoles 
have been discussed in connection with 'topological inflation' [^, ^. More recently per- 
turbations of static global monopoles were studied in the context of 'critical behaviour' 
in spherically symmetric gravitational collapse 0. 

In the present paper, we shall derive some of the qualitative and quantitative results 
of from the field equations - a set of non-linear ODEs - obeyed by these solutions. As 
described in there is a one-parameter family of static global monopoles parametrized 
by the Higgs vacuum value r] (taken in units of Mpi = 1/a/G) for < 77 < ?7max ~ ^J3/8'K. 

As long as ?7 < ^l/Sn the space-time of these solutions has a deficit solid angle at infinity 

becoming equal to Att for rj = ^1/8tt. Solutions with y^l/Svr < r] < rj^^x were found 
to have a cosmological horizon, outside of which the solutions oscillate with decreasing 
amplitude about their asymptotic values. As will be demonstrated, this behaviour can 
be readily understood from the field equations. 

Our main result is a simple derivation of the value r^max as being exactly equal to 

^S/Sn. At this value of t] the numerically determined global monopoles are found to 
bifurcate with the de Sitter solution obtained for vanishing Higgs field - the non- vanishing 
vacuum energy yielding a cosmological constant. 

A linear stability analysis of both types of solutions shows the stability of the global 
monopole for all values < rj < r^max, whereas the de Sitter solution (considered as a 
solution of the Einstei n-Higg s system) is stable against perturbations with support inside 
the horizon for r] > ^JsJstt and unstable for smaller values of r] with an accumulation 
of unstable modes for 77 — > 0. The stability of the global monopole is in agreement with 
numerical studies of time-dependent solutions performed in 

Following the notation of we put (j)"" = f{r)f"' for the Higgs field and 

ds^ = -A^fidt^ + — + r^dn^ (1) 
/i 

for the spherically symmetric line element in Schwarzschild coordinates. The resulting 
static field equations are (we prefer to use rationalized units putting / = \/47rf, fj = 
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Solutions with a regular origin obey the boundary conditions 

f{r) = ar + 0{r^), V^(r) = a + ©(r^), /i(r) = 1 + ©(r^), A(r) = + ©(r^) (6) 

and are uniquely specified by the choice of a = ip^O) and Aq. 

Such solutions stay regular as long as /x > 0. Their global behaviour is characterized 
by three possibilities: 

1. /i has a zero for finite r = ro with diverging ip, while A, / and ^/JI1p stay finite. This 
type of singularity is just a coordinate singularity related to the stationarity of r 
considered as a metrical function. In terms of a 'geodesic' coordinate r = Iq 
the function r(r) reaches a maximum Tq = r(ro) and then runs back to r = 
developing a curvature singularity there. This type of behaviour is generic for 
solutions with a regular origin. 

2. fi has a zero for finite r = rh with finite ^Z', / and A. Such solutions possess a 
cosmological horizon at r = rh enforcing the boundary conditions 
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with /h = /(rh) etc.. 
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3. ;U > for all r > 0, i.e. solutions staying finite for all r. The asymptotic behaviour 
for r — > oo is given by 

(9) 



/i(r) = 1 - 2f ^ - + O(^) 
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and some constants c and d depending on the solution. 
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The asymptotic behaviour of /i entails a 'conical' singularity at r = oo correspond- 
ing to deficit solid angle A = Snff. For ff = \ this deficit angle becomes An and 
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the character of the solution for r — > oo changes to 

Mr) = ^(l + 0(i)) (11) 
^"<'>= = 71-7? + 0(i)+c.— "(l + 0(i)) (12) 



with 1^ = ly I and again constants c and d depending on the solution. 

For fj > 1/a/2 the asymptotic behavior for r — > oo becomes oscillatory due to the 
sign change of a;^. However, the factor 1/r in front of the exponential leads to a 
damping of the amplitude and /(r) — s> f/ for r ^ oo as is clearly visible from the 
numerical results of [0]. 

Both the solutions of type 2) and 3) require fine-tuning of the parameter a char- 
acterizing solutions regular at r = 0. The numerical results of indicate that for 
< fj < f/max exactly one such solution with non-vanishing Higgs field exists. The cor- 
responding graph of a{f]) is shown in Fig.|^. In the following we shall try to explain the 
essential features of a{f]). 

For small values of f] the self-gravitating monopole can be understood as a pertur- 
bation of the flat monopole, which is obtained for 77 — > after a rescaling f = f]f and 
r = f/fj solving the equation 
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(rY) =/[2 + -(/^-l)J. (13) 

Solving this equation numerically with the boundary condition /(f) = af + O(f^) one 
finds a globally regular solution for a ~ 0, 3578. Expressed in terms of / we get 

f{r)=ffdr + 0{r^) = ar + 0{r'^) (14) 

and thus a = ffa. This behavior is clearly visible in FigJ^. 

We also see from this figure that the function a{ff) after reaching a maximum turns 
back to a = with diverging derivative. The latter may be understood as a consequence 
of the invariance f —f of the field equations in combination with the smooth behaviour 
of a{fj) at r/max- Next we will use this "empirical" observation to derive the exact value 
for r^max- 

Suppose we take a globally regular solution f{a{f]),f],r) for some fixed value of fj 
and vary a{fi) — > a(f/) + 6a. Then generically we will get a solution of type 1). However, 
if da{ff)/df] = 00 an infinitesimal change of a will carry us to another regular solution 
for the same value of fj. Thus df /da will be bounded for r 00 resp. r — rh for points 
where da/dfj = 00 in contrast to points where this derivative is finite. As seen from 
Fig.H the solution for fj = f^max corresponds to a(f/) = and thus /(//max, r) = 0. This is 
nothing but the de Sitter solution given by integrating Eq.(^) 

fi'dsir) = 1 - (— ) with Th = . (15) 

Figure |^ shows the approach of the monopole solutions to the de Sitter solution for 
V Vma.x- According to what was said above, the de Sitter solution must have the 
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Figure 2: a = /'(O) for the static monopole as a function of ff. Observe that a — ^ 
as f/^ — i> 3/2 signifying the approach to de Sitter space. Note that the critical point in 
a does not occur ai ff = 1, which was found to be the critical point in the core radius 
in i. 
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Figure 3: Demonstration of approach to de Sitter solution. The field /u'(r) is shown 
as ff approaches 3/2. The solid line denotes the de Sitter solution for r^^ = 3/2 while 
the other profiles show jj! for the static monopoles. The solutions approach the de Sitter 
solution for an increasing domain of r until ultimately driven to oscillate around the 
regular solution /i' = 0. 
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bounded "zero mode" df /da for rj = ?7max- Linearizing the field Eqs.(0-|^) on the de 
Sitter background we obtain (putting ip = r6f) 



^4 ^2 



W)'-((l ->')'= (16) 



Introducing x = r/r^ we get 



The solution with the correct boundary condition = + O(x^) for a; ^ valid for 
if = df /dais given hy ip = x^ obtained for f/^ = 3/2. 



Thus we have found a simple derivation of r^^ax = Since the r.h.s. of Eq. (p!7D 



becomes more negative for ff < 3/2 the solution (f{x) fulfilling the correct boundary 
condition at x = develops a zero before reaching the horizon at x = 1. According 
to the "Jacobi criterion" this implies an instability of the de Sitter background for 
f] < Linearizing the time-dependent field equations around the static solution we 
let 

/(r,t) = /o(r)+e^-*5/, (18) 

where /o represents the static Higgs field. Unstable modes of the static solution corre- 
spond to modes with a;^ < 0. The Jacobi criterion relates the existence of such negative 
modes to the existence of zeros of the solution for uj^ = (zero mode) with the correct 
boundary condition at r = 0. 



For V > y ^ the solution (p{x) remains positive for all x indicating the stability of 
the solution under spherically symmetric perturbations with support inside the horizon. 
On the other hand for f/ — ^ the de Sitter solution develops more and more unstable 
modes manifested by the existence of bounded zero-modes with more and more zeros. 
The latter are readily obtained with the polynomial ansatz 

^K{x) = j2ckx'''- (19) 
fc=l 

Inserting this ansatz into the Eq.(p!7|) yields the recursion relation 

2k{2k + l)-2-^ 
= {2k + 2){2k+l)-2 

and the condition (derived from ck+i = 0) 

3 

= Ki2K + l)-l ^^^^ 

yielding f]f = 3/2, f]^ = 1/3, fj^ = 3/20 etc.. 

Eq.(^) can be transformed to the hypergeometric differential equation as follows. 
Putting ^ = Yzfj and (p{x) = ^{1 + with 2m{2m — 1) = 3 + 6/r/^ we obtain 

+ 1)/ + ((2m + 3)e +l)y' + {m + ify = (22) 



6 



4 



2 



2 



4 










1 


.51 - 


7]^ = 


1 


.50 J 


T]^ = 


1 


.49 - 


77^ = 


1 


/3 - 


T]^ = 


3 


/20^ 



40 6( 



Figure 4: Perturbations to the de Sitter solution. For ff > 1.5, the perturbation is ev- 
erywhere positive indicating linear stability. For fj'^ < 1.5, zero-crossings exist indicating 
instability of the de Sitter solution. Additional zeros appear as f/^ passes through the 
values fj]^ found from Eg. (pT]) , the first few of which are shown here. 
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Figure 5: Perturbations to the static global monopole. For all r]^, the perturbation is 
everywhere positive indicating linear stability. 



fulfilled by the hypergeometric function F{m + 1, m + 1, 5/2, ^) P], the polynomial solu- 
tions (Pk{x) corresponding to the Jacobi polynomials yxiO- 

The stability of the fiat global monopole (under spherically symmetric perturba- 
tions) implies its stability for small values of f/. According to standard arguments for 
one-parameter families of solutions we expect no change of stability up to the bifur- 
cation point f] = f/max with the de Sitter solution. In order to analyze this hypothesis 
we have performed a detailed numerical stability analysis for the global monopole using 
again the Jacobi criterion. For solutions with a horizon our stability analysis is restricted 
to the region inside the horizon. This is justified by the fact that no perturbations can 
penetrate through the horizon from the region where /i < (similar to the case of 
stability analysis performed for black holes). 

As is well-known ||^ the metric perturbations can be integrated out and one obtains 
a Schrodinger type perturbation equation for ip. Introducing the coordinate a through 
^ = -fo one gets 

-^ + V^ = (23) 
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with 



r2 1 
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(24) 



and the boundary condition ^p{r) = cr^ + 0((T^) for cr — 0. Our numerical analysis 
shows (see Fig.^ that the solutions ip have no zero for all < f/ < ^3/2 confirming our 
hypothesis about the stability of the global monopole. 



1 Acknowledgments 

This work began at the ITP at the University of California, Santa Barbara on occasion of 
the workshop on Classical and Quantum Physics of Strong Gravitational Fields and has 
been supported in part by the National Science Foundation under Grant No. PHY94- 
07194. The hospitality of the ITP is gratefully acknowledged. SL is also thankful for the 
financial support of Southampton College. 



References 

[1] S.L. Liebling, "Static gravitational global monopoles" , 
t^r-qc/9906014| . 

[2] A. Linde, Phys. Lett. B 327 (1994) 208. 

[3] A. Vilenkin, Phys. Rev. Lett. 72 (1994) 3137. 

[4] S.L. Liebling, Phys. Rev. D 60 061502 (1999). 

[5] I.M. Gelfand and S.V. Fomin, Calculus of Variations, 
Englewood Cliffs, N.J., Prentice Hall, 1963. 

[6] A. Erdelyi, M. Magnus, F. Oberhettinger and F.G. Tricomi, Higher Transcendental 
Functions, Vol.2, 
McGraw Hill, New York, 1953. 

[7] K. Maeda, T. Tachizawa, T. Torii and T. Maki, Phys. Rev. Lett. 72 (1994) 450. 

[8] H. Hollmann, Phys. Lett. B 338 (1994) 181. 



9 



